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ABSTRACT 

The discovery of two neutron stars with gravitational masses ~ 2 Mq has placed a 
strong lower limit on the maximum mass of nonrotating neutron stars, and with it a 
strong constraint on the properties of cold matter beyond nuclear density. Current upper 
mass limits are much looser. Here we note that, if most short gamma-ray bursts are 
produced by the coalescence of two neutron stars, and if the merger remnant collapses 
quickly, then the upper mass limit is constrained tightly. If the rotation of the merger 
remnant is limited only by mass-shedding (which seems probable based on numerical 
studies), then the maximum gravitational mass of a nonrotating neutron star is ~ 

2 — 2.2 Mq if the masses of neutron stars that coalesce to produce gamma-ray bursts 
are in the range seen in Galactic double neutron star systems. These limits would be 
increased by ~ 4% in the probably unrealistic case that the remnants rotate at ~ 30% 
below mass-shedding, and by ~ 15% in the extreme case that the remnants do not rotate 
at all. Future coincident detection of short gamma-ray bursts with gravitational waves 
will strengthen these arguments because they will produce tight bounds on the masses 
of the components for individual events. If these limits are accurate then a reasonable 
fraction of double neutron star mergers might not produce gamma-ray bursts. In that 
case, or in the case that many short bursts are produced instead by the mergers of 
neutron stars with black holes, the implied rate of gravitational wave detections will be 
increased. 

Subject headings: dense matter — equation of state — gamma rays: bursts — gravita¬ 
tional waves — stars: neutron 


1. INTRODUCTION 


The state of cold matter beyond nuclear density cannot be resolved strictly with laboratory 
experiments, and nuclear theories diverge strongly in their predictions for such matter. Thus astro¬ 
nomical observations are sought for guidance. Important constraints were obtained from the dis¬ 
covery of stars with gravitational masses of M = 1.97±0.04 Mq (PSR J1614-2230; iDemorest et al 
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2010 ) and M = 2.01 ± 0.04 M 0 (PSR J0348+0432; Antoniadis et al. 20131 ). Radius 


measurements 

would be helpful, but current estimates are dominated by sys tematic errors although hope exists 


for future X-ray and gravitational wave measurements (see [Miller] 2013 for an extensive discussion) 


It would also be useful to have an upper limit to the maximum gravitational masses of nonrotating 
stars, in addition to the_current ) > 2 M 0 lower limit (note that PSR J0348+0432 has a spin 

period of 39 ms Antoniadis et ah 2013], which is long enough compared to the ~ 1 ms minimum 


Rhoades & RufRni 

(1974 

) found Mmav 3.2 M m , and 

Kalogera & Bavm 


(199Q) performed an updated treatment that found M max <j 2.9 M 0 , where the differences depend 


primarily on the density up to which we believe we know the equation of state (EOS) of cold matter. 

Ongoing observations of short gamma-ray bursts, and recent theoretical considerations of their 
mechanism, might provide a route to tighter upper limits on the maximum mass. Short gamma-ray 
bursts have long been thought to be most pro bably caused by the merger of either two neutron 
stars or a neutron star and a black hole (see IBer ger 2014 for a recent review). From the stand¬ 
points of energetics and time scales other candidates exist, such as the accr etion-induced col lapse 


of a white dwarf followed by magnetar-like rapid spindown (e.g., iDar et al.lll992l : lYi &; Blackman 


19971 : iMetzger et al.l 120081 ) . However, in additio n to the burst r a te being roughly consistent with 


the expected compact object merger rate (e.g., Fon g et al 


2012 


Wanderman &; Piran 


201 


5), and 


magnetar models being disfavored by the lack of confirmed radio transients (Metz ger et al. 120151 ). 
the merger picture has received rec ent support from estimates of the spat i al offsets of bursts 


from their most like l y galactic hosts dBelczvnski et al.ll200fil : iFong et al.l 120101 : IChurch et al.l 12011 


Fong &; Bergeri l2013: iB ehro ozi et al. 2014])- These offsets are broadly consistent with expectations 
based on the supernova recoil that accompanies the births of neutron stars and black holes, and 
on the time needed after formation for the binary to coalesce due to the emission of gravitational 
waves. 

The compact object coalescence model for short gamma-ray bursts has been explored numer¬ 
ically and analytically with progressively greater fidelity over the last several years. Observations 
of afterglow spectral changes characteristic of jet brea ks imply that short bursts, like their long 
counterparts, have large bulk Lorentz factors T > 100 (IKrolik fc Pierl Il99ll : iFenimore et al.lll993l : 
Baring fc Harding! 1997). Thus the jets that we see from short bursts must avoid being loaded with 
too many baryons, because this would slow down the jet (ISherni & Piran 1990) and would delay 
and lengthen the burst. If the compact objects are a black hole and a neutron star, it is believed 
that this happens automatically because the orbital axis will naturally be nearly free of baryons 
due to the existence, from the beginning, of an event horizon in the system. 

However, if the com pact o bjects are both neutron stars then contamination of the jet by 
baryons is possible. iMurguia-Berthier et al. (2014j) argue that in order to produce the observed 
high Lorentz factors, the merger remnant must collapse within no more than 100 milliseconds 
after the initial merger, otherwise baryons driven outwards by interactions with neutrinos would 
overload the jet and cause the burst to last longer than is typical for short gamma-ray bursts (see 
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also Rosswog Sz Ramirez-Ruiz 2002). In this sense, it is similar to the argument that the duration 
of a long gamma-ray burst is comparable to or larger than the time needed for the jet to escape 
the envelope of the massiv e star who se co ll aps e caused the burst. For an alternate view see the 
“time-reversal” scenario of Ciolfi & Siegel 2015, who suggest that delayed collapse of a uniformly 
rotating star can explain both the prompt gamma -ray e mission and the occasional 10 2 — 10 5 s 
X-ray emission that follows; note, however, that Margalit et ah 2015 argue that the collapse of a 
uniformly rotating star will not produce a disk and therefore will not generate a powerful jet. More 
work is clearly required. 


If this picture is correct and if most short gamma-ray bursts indeed come from mergers of two 
neutron stars rather than from the merger of a neutron star and a black hole, then we can infer an 
upper limit to the masses of nonrotating neutron stars. The essence of the idea, which we elaborate 
in Section [2l is that EOSs that predict large maximum masses for nonrotating stars also predict 
that the combined mass of two neutron stars can be supported stably in a uniformly rotating con¬ 
figuration, which would mean that short gamma-ray bursts are not producecJE Assuming that the 
double neutron star systems we see in our Galaxy are representative of the population that pro¬ 
duces short gamma ray bursts, we show that this leads to tight constraints on the maximum mass. 
We find that if the uniform rotation of merged remnants is limited only by mass-shedding (which is 
consistent with published neutron star merger simulations) then the maximum gravitational mass 
for nonrotating stars is ~ 2 — 2.2 Mq, depending on the masses of the double neutron stars that 
produce gamma-ray bursts and somewhat on the EOS class that we consider. Although remnant 
rotation at the mass-shedding limit is consistent with all current simulations, we also consider 
more slowly rotating remnants, which are conceivable if a large fraction of the angular momentum 
is removed from the system very soon after merger. In the extreme case that the remnant has no 
angular momentum our upper mass limits are increased by ~ 15%. 


These tight limits would have important implications for cold matter beyond nuclear density. 
Depending on the range of masses of neutron stars in mergers, this could also imply that a fair 
fraction of mergers do not lead to short gamma-ray bursts. This raises the question of what form 
their still-large energy release would take. It also implies that the merger rates inferred from short 
gamma-ray bursts need to be increased. This would increase the likely merger detection rates for 
ground-based gravitational wave detectors such as Advanced LIGO , Advanced Virgo, and KAGRA 
( Doolev et al. 2014 ). 


If instead short gamma-ray bursts are produced by the coalescence of a neutron star with a 
black hole, then our argument does not apply and we therefore cannot use this argument to place 
an upper limit on the maximum mass of a nonrotating neutron star. This is because, in that case, 
an event horizon already exists and the jet funnel will therefore presumably already be clean. The 
news would then be even better for ground-based gravitational wave detectors, because given that 


x As we were performing the work for this paper we learned that a s imilar idea for lim iting neutron star maximum 

(2Qll)- 


masses, with somewhat different methodology, was being prepared by Ffyer et .al 






















-4- 


black holes are more massive than neutron stars, the signal produced will be stronger and thus 
visible to greater distances than the signal from a double neutron star merger. Therefore, for a 
given rate of short gamma-ray bursts per volume, the detection rate will be significantly greater 
than for double neutron star mergers. 

In Section [2] we discuss our assumptions and method of analysis. In Section [3] we give our 
results. We motivate the parametrized forms of our EOSs and the allowed ranges for those param¬ 
eters, and then give the mass limits that result. We finish in Section []] by discussing how future 
electromagnetic observations of gamma-ray bursts that are also detected with gravitational wave 
instrumentation can make our results more rigorous. In particular, we show that even if only the 
chirp mass, rather than both masses independently, can be measured, the uncertainty about the 
separate masses has only a small effect on our constraints. 


2. METHODS 


Our pri mary tool for analy sis is the publicly available Rotating Neutron Star (rns) code 
developed bv IStergioulas fc Friedman! (1995). This is a flexible code that computes the structure 
and external spacetime of an axisymmetric, uniformly rotating star for a given EOS. We have 
modified the code slightly so that for a given EOS it can output (1) the maximum gravitational 
mass for a nonrotating star, (2) the baryonic rest mass corresponding to a given gravitational 
mass for a nonrotating star, and (3) the maximum baryonic rest mass for a uniformly rotating 
configuration, which is limited by either mass-shedding alone or, in addition, by a limit on the ratio 
T/W of the rotational kinetic energy to the gravitational binding energy (note that T/W at mass 
shedding typically ranges from ~ 0.1 for very soft EOS to ~ 0.15 for very stiff EOS; see, e.g., Table 5 
of Cook et ah| 1994). One way in which a threshold on T/W could be relevant is if the merged 
remnant develops a nonaxisymmetric instability, because in principle such nonaxisymmetries could 
lead to the emission of gravitational waves that would reduce T/W to the threshold value. However, 
treatments of relativistic fluids suggest that the threshold value js T/W > 0 .2 for realistic neutron 
star compactnesses (Shapiro &; Zane 1998; Gondek-Rosinska & Gourgoulhon 2002), which is well 
above the mass-shedding limit. Thus secular instabilities are not likely to be important, so we 
consider T/W thresholds only as a way to determine the effect on our limits of, e.g., rapid angular 
momentum loss due to magnetic braking. 

We have tested our implementati on of the code by comparing our an swers to (1) — (3) a bove 
with those given in Co ok et al.J (1994) for several tabulated EOSs and in Read et ah (20091) for 
several parametrized approximations to other EOSs. We find agree ment to better than about 1% 
in all cases except for the very soft EOSs F and G used in ICook et ah) (1994|); these EOSs are ruled 
out by the existence of neutron stars with gravitational masses ~ 2 Mq . We note that for most of 
the EOSs listed in Cook et ah (1994), the sound speed becomes superluminal (c 2 = dP/dp > c 2 ) 
above some energy density that is usually less than the central energy density for the maximum 
mass rotating and nonrotating configurations. If we impose a c s = c upper limit then the maximum 
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masses often decrease significantly: for example, for EOS UU, the maximum gravitational mass for 
a nonrotating star drops from M max = 2.20 M@ to 2.06 Mq and the maximum baryonic rest mass 
for a uniformly rotating star drops from Mb a ry,rot,max = 3.12 Mq to 2.94 Mq. 

Our limits are conservative in the sense that we assume that the remnant is cold, and that 
it rotates uniformly rather than differentially. A hot remnant can support more mass than a cold 


remnant (see, e.g., the recent treatment in K aplan et al . 2014L and_differential rotation can sup 


port more mass than uniform rotation (e.g., Baumgarte et ah 12000 ). However, we expect that 


temperatures comparable to the Fermi temperature ~ 10 12 K, which are required to provide sig¬ 
nificant extra support, will exist for at most an extremely short tim e becau se neutrino emission 


will remove the energy efficiently. We also expect, as suggested by Shapiro (2000), that internal 


magnetic fields will be amplified rapidly enough by differential motion that the angular momentum 
will be redistributed into a state of uniform rotation. If either of these assumptions is incorrect 
then the upper limit to neutron star masses will be lowered somewhat. Our upper limit would also 
be lowered if some of the baryonic mass in the two stars ends up in a disk or outflow rather than 


amounts to on 

y few x 10 2 Mq (e.g., 

Shibata & Urvu 

-— 

2000: Shibata et al. 

2003: Duez et al. 

2004 

Liu et al. 2008 

Kiuchi et al. 2009; Rezzolla et al. 

2010 

Fan & Wei 

2011; 

Giacomazzo et al. 

2013 

Bernuzzi et al. 

2014; Kastaun & Galeazzi 2015; Dietrich et al. 2015 

)• 


Our assumptions might not be conservative if the dynamics of neutron star mergers are such 
that the remnan t rotates at a rate less than the mass-shedding limit for the remn ant. Current 
simulations (e.g.. iBernuzzi et alJl2014l : iKastaun fc Galeazzill2015l : iDietrich et alJl2015l ) suggest that 
the total angular momentum of the remnant plus disk and outflows does exceed the mass-shedding 
limit, so we consider this our standard case. However, if the rotation limit is tighter, less mass can 
be supported by the spinning remnant than in our assumptions, so a larger range of EOSs would 
satisfy our criterion for short gamma-ray bursts, and hence the maximum allowed gravitational 
mass for a nonrotating star would be larger than the mass we derive in our standard analysis. 
Our mass limit would also be increased if the remnant has a strong enough poloidal magnetic field 
(either initially, or developed via, e.g., a dynamo induced by differential r otat ion) that it _spins 
down significantly during the ~ 0.1 s allowed by the argument of Murguia-Berthier et ah (2014). 
This would require a poloidal field with an extremely large characteristic surface strength of at 
least B ~ few x 10 16 G. To take these possibilities into account we also explore the mass limits that 
come from the assumption that (as an arbitrary, round number) T/W = 0.1 limits the rotation, 
and the mass limits that would apply in the extreme case that the remnant is nonrotating. 

Our method is the following: 


1. Start with an assumed EOS and two neutron stars, which we assume to be nonrotating 
(because the f astest-sp inning neut ron star in a double neutron star system has a frequency 
of only 44 Hz Burgav et al.l 1200311 and tid a l torques will no t spin the stars up significantly 


even near merger 


Bildsten &; Cutled Il992l : iKochanekl [1992] ). Let the gravitational masses 
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be Mi and M 2 , and let their corresponding baryonic rest masses for the chosen EOS be 
Mbary.i and Mb aryi 2 - We investigate three pairs of masses: Mi = 1.25 M 0 and M 2 = 
1.35 M 0 (compara ble to PSR J0737-3039A,B, which is the the lightest double neutron star 
pair yet discovered Burgay et al. 2003; Lyne et ah 2004 ]): Mi = M 2 = 1.35 M 0 (comparable 
to the average gravitational mass of double_ neutron stars discovered in our Galaxy, and 


simil ar to the systems PSR B1534+12 Stairs et al.ll2002l] and PSR B2127+11C Jacobv et al 


2006]); and Mi = 1.35 M 0 and M 2 = 1.45 M 0 (c omparable t o PSR B1913+16. wh ich is 


the heaviest double n eutron star pair yet discovered [Tavloiill992l : IWeisberg et al.ll2010l ]h See 


Kiziltan et al. (2013) for a recent summary of the masses in double neutron star systems. 


2. If, when they merge, the stars produce a gamma-ray burst, then by the logic of Murguia-Berthier et al 
(2014j) the remnant must collapse quickly to form a black hole. Thus the baryonic rest mass 
M rern bary of the remnant must exceed the stable limit of a uniformly rotating neutron star. 


3. We assume that M rem> b ary = Mb aryjl + Mmry, 2 > which is the maximum possible. Any matter 
that goes into an outflow or a disk that lasts for more than 0.1 s will reduce the remnant 
mass and strengthen our argument. 


4. We compare M remi bary with the maximum baryonic rest mass Mb ary ,rot,max that can be sus¬ 
tained by a uniformly rotating star for the assumed EOS. The rotation is limited by either 
(a) mass-shedding, or (b) a limit on T/W, whichever is more restrictive. Our primary re¬ 
sults are based on (a), which we believe to be the most realistic case, but we also explore 
T/W = 0.1 and T/W = 0 (nonrotating). 

5. If Mbary,rot,max < M rem) b a ry then the remnant collapses and the equation of state is viable. 
For this equation of state we can also compute the maximum gravitational mass M max for a 
nonrotating star. 

6. We therefore search the parameter space in our EOSs to find the largest M max that is viable 
by our short gamma-ray burst criterion. This is the number we report. 


3. RESULTS 

3.1. Parametrized equations of state 

The neutron Fermi momentum in pure neutron matter at densities below twice nuclear sat¬ 
uration density (n < 2 n s , where n s = 1.6 x 10 38 cm -3 ; the corresponding mass density is p s = 
m. n n s = 2.7 x 10 14 g cm' 3 ) is less than ~ 420 MeV/c. At n = 2 n s , two neutrons with opposite 
momenta on the top of the Fermi sphere have a total energy of about 320 MeV (in the lab frame) 
or 160 MeV (in the center of mass frame). This is barely above the pion production threshold of 
~ 140 MeV and indicates that a treatment of dense neutron matter based on neutron degrees of 
freedom interacting through a potential should be adequate. The character of interactions between 


























- 7- 


the nucleons can be inferred from elastic nucleon-nucleon scattering data and the spectroscopy of 
light nuclei {A < 20). Modern many-body methods, whether computational or analytical, can then 
be used to infer the zero temperature EOS. This program has been carried out in the last few years 
and represents an important step towards an understanding of dense matter from first principles, 
at least at the low end of densities relevant to neutron stars. 


There are two versions in th e literature of the program sketched above. The first one (jHebeler et al 


2010l : iHebeler &; Schwenkl 1201(1) attempts to describe the nucleon-nucleon interaction using effec¬ 
tive field theory ideas ( Weinberg 1991 : Bedaaue fe van Kolck 2002 : Beane et al. 2000 : Phillips 2002 : 
Epelbaum et ah 20091) to extract the nucleon-nucleon potential from quantum chromodynamics. In 
this approach there is a systematic expansion of the interaction in powers of the momentum and 
pion masses and uncertainties can be quantified a priori. Three-body forces appear as a small but 
important effect. Unfortunately, this approach describes the nucleon-nucleon phase shifts up to 
relatively small momentum and is adequate only for densities below nuclear saturation. 


A second approach is more useful for our purposes (Gandolfi et ah 2012, 2014), and it forms the 
basis for our primary parametrized equation of state, which we call EOS1. In this approach, phe¬ 
nomenological nucleon-nucleon potentials are fit to nucleon-nucleon elastic scattering data. These 
potentials describe the phase shifts well up to energies around 600 MeV (in the lab frame). Two- 
nucleon potentials are, however, not enough to describe matter even below nuclear saturation 
density. Effective theory arguments as well as studies of light nuclei demonstrate that three-body 
forces are required. The importance of the three-body forces increases with density and is sub¬ 
stantial at n = 2 n s . Some components of the three-body force can be extracted by fitting light 
nucleus energy levels. However, in neutron matter, unlike in nearly symmetric matter, only the 
isospin-3/2 channel is relevant, so the three-body components that are derived from light nuclei are 
not the dominant ones in neutron matter. Instead, Gando lfi et al. (2012, 2014) used a variety of 
three-body forces with differing functional forms and ranges varying over a factor two around the 
pion Compton wavelength, where the strength of each fo rce w as fixed so that the symmetry energy 
lies within the empirically observed range (|Tsang et al.ll2012l ). These interacti ons were used to ob¬ 
tain the EOS using the auxiliary field Green’s function Monte Carlo method Schmidt fe Fantoni 
1999). The errors arising from the Monte Carlo and infinite volume extrapolation are negligible for 
our purposes. The three-body force is, then, the largest source of uncertainty about the neutron 
matter EOS at densities below 2 n s . 


The results of Gan dolfi et alj (2012, 2014) are well fit with the convenient parametrization 


e(n) = n m n + a(n/n s ) a + b(n/n s y 


(1) 


for the mass-energy density e as a function of the number density n. The parameters a and a 
depend primarily on the two-body force and are well constrained by scattering data to lie in the 
ranges 12.6 MeV < a < 13.0 MeV and 0.47 < a < 0 .50. The parameters b and /3 are more sensitive 
to the three-body force; in the analysis of Gandolfi et al. (2012, 2014]) they vary over the ranges 





































































3.2 MeV < b < 5.2 MeV and 2.1 < (3 < 2.5 q However, in order to be conservative, we will present 
results obtained by doubling the uncertainty in these parameters. Thus the ranges we search are 
12.4 MeV < a < 13.2 MeV, 0.45 < a < 0.52, 2.2 MeV < b < 6.2 MeV, and 1.9 < 0 < 2.7. 


For EOS1 we adopt the equation of state given above to a threshold rest mass density pthresh.i, 
where we explore the range 1.7 p s < pthresh.i <2 p s . Between /?thresh,i and pthresh ,2 (which we do 
not constrain except to require that it exceed pthresh.i) we assume a constant sound speed c s that we 
allow to be anywhere between c s = c/2 and c s = c (the low-density equation of state in Equation |T] 
is itself never acausal). We assume that c s = c above /9thresh,2 because the central mass-energy 
density of a maximum mass nonrotating star is greater than that of a maximum mass uniformly 


r otat ing star (compare, e.g., the central densities of the nonrotating stars i n T able 4 of Cook et al 


1994] with the central densities of the rotating stars in Table 5 of Cook et al. 19941) • Therefore a 
transition above some density to the hardest possible EOS can increase the maximum mass of a 
nonrotating star without affecting the maximum mass of a uniformly rotating star. 

The class of equations of state included in EOS1 is very large and was designed with the 
objective of maximizing the maximum non-rotating mass while allowing for short gamma-ray bursts 
to follow mergers. Still, it does not include some possibilities that, while not favored, are not 
excluded by hard evidence. The calculations leading to the low density part of the EOS1 assume 
that nucleons (and electrons and muons) are the only relevant deg rees o f freedom at densities below 
1.7n s . This expectation can be frustrated if a pion condensate (jMigdall 1978!), or hyperons, become 
important at these low densities. Equations of state with any extra degrees of freedom besides 
nucleons tend to be too soft to support a maximum non-rotating mass of 2 Mq but the possibility 
remains that the correct equation of state is softer than EOS1 at n < n s but stiffens quickly at 
higher densities. A similar situation would also be obtained in th e ev en more unlik ely poss ibility 
that a transition to quark matter occurs at densities below 1.7 n s ([Madsen 1999; Alfor d et ah 20081). 
Finally, a very different scenario arises if strange quark matter is the true ground state of matter at 
arbitrarily low de nsities. In this case the low density equation of state would have a very different 
form from EOS1 ([Bodmer 1971; Witten 1984) and our results would not apply at all. 


Our secondary parametrized_EOS, which we call EOS2, is a slight modification of the piecewise 
polytrope introduced by [Read et al. (120091 ). EOS2 represents the pressure as a function of density 


2 There is a strong correlation between the parameters b and /3 that further constrains the set of equations of state. 
In the spirit of being as conservative as possible we will neglect this correlation in our work. 
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with four parameters po, Ti, T 2 , and Ta: 


P 

oc p 1 ' 58425 

p < 2.44034 x 10' g cm 3 


oc p 1 - 28733 

2.44034 x 10 7 g cm -3 < p < 3.78358 x 10 11 


oc p 062223 

3.78358 x 10 11 g cm" 3 < p < 2.6278 x 10 12 


= 3.594 x 10 13 dyn cm -2 p 1 - 35692 

2.6278 x 10 12 g cm -3 < p < po 


oc p Fl 

po < p < 10 14 ' 7 g cm -3 


oc p V2 

10 14 ' 7 g cm -3 < p < 10 15 g cm” 3 


oc p Fa 

10 15 g cm” 3 < p . 

The 

normalization in the 2.6278 x lO 1 ^ 

g cm” 3 < p < po range is taken from 


a ” 3 

3 


( 2 ) 


Read et al. (j2009 l: note that their table lists the normalization divided by 


We impose the additional limitation that if there is any density p caus at which these expressions 
would pre dict that the EOS becomes acausal above some density / 0 caus , we set c s = c at p > p caus (see 
Ellis et al .1 2007 for a discussion of the causality limit). That is, we do not constrain the parameter 
space (po, Ti, T 2 ,1^) based on the requirement of causality; instead, if a given combination predicts 
c s > c above some density, we set the sound spe ed to be e qual to c above that density. This 
approach generalizes somewhat the EOS of Read et al. (2009), and is consistent with 
our philosophy of considering as broad a set of EOSs as possible. The prefactors for the 
pressure in each density range are set by pressure continuity at each of the density boundaries, 
which are anchored by the pressure in the 2.6278 x 10 12 g cm' 3 < p < po range, and all the 
densities are measured in g cm -3 . The ranges that we search in the parameters are po = (1 — 8) x 
10- 4 M q /(GM q /c 2 ) 3 = (1 - 8) x 6.173 x 10 13 g cm' 3 ; T x = 1.5 - 5; T 2 = 1 - 5; and T 3 = 1 - 5. 


3.2. Mass limits 


We use Powell’s direction set method (jPowelll 1964) to maximize the maximum gravitational 
mass of a nonrotating star over the parameters in each of our EOS models. For both EOSs, 
the maximum mass is nearly unimodal over the parameter space, which makes searches relatively 
smooth and reproducible. 


We display our results in Figure [D For each of the three combinations of gravitational mass 
(1.25 M q , 1.35 Mg), (1.35 Mg, 1.35 Mg), and (1.35 Mg, 1.45 Mg) we show the mass limits we 
obtain for EOS1 and EOS2 and for different rotation limits: mass shedding only (no T/W limit); 
the more restrictive of mass shedding and T/W < 0.1; and a nonrotating remnant. This figure 
shows that if short gamma-ray bursts are produced by mergers of double neutron stars similar to 
those we have discovered in the Galaxy and if the remnant is unstable in a uniformly rotating 
state, then the limits on the maximum mass of a nonrotating neutron star are extremely tight, 
particularly if mass shedding alone sets the limit on uniformly rotating merged remnants. We see 
that the limits from EOS1 are tighter than the limits from EOS2, because EOS2 allows greater 
freedom for the EOS at high densities. We also see that if low-mass double neutron star mergers 
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produce short GRBs then the maximum gravitational mass of a nonrotating neutron star is already 
known with high precision. 


4. IMPLICATIONS AND SUMMARY 

We have shown that if short gamma-ray bursts are produced when mergers of two neutron 
stars lead to a rapid collapse to a black hole, and if the neutron star masses are similar to what 
we see in Galactic double neutron star systems, then the maximum mass of nonrotating neutron 
stars is constrained tightly from above as well as from below. This is because if the maximum mass 
exceeds our limit, the merger of two neutron stars of the specified masses will produce an object 
that remains stable even after internal angular momentum transport produces uniform rotation. 
Thus prompt collapse would not happen; any collapse would be delayed by at least as long as it 
would take to radiate angular momentum from the system. 

Much more certain limits on the maximum mass will be obtained when ground-based grav¬ 
itational wave detectors see the gravitational radiation from double neutron star coalescences in 
coincidence with short gamma-ray bursts. This is because such coincidences will allow us to limit 
strongly the total baryonic rest mass of the merger. The smallest such total mass associated with a 
successful burst will place the tightest upper limit on the maximum mass of a nonrotating neutron 
star. 


To determine how tightly we could constrain the total baryonic rest mass in a merger, we note 
that the mass-related quantity that will be most precisely measured from gravitational wave trains 
is the chirp mass 


M ch = rf' b M tot 


(3) 


where Mt 0 t = + M 2 is the total gravitational mass and q = M\M 2 /M^ ot = q/( 1 + q ) 2 is 

the symmetric mass ratio, where q = M 1 /M 2 < 1 is the standard mass ratio. The chirp mass is 
essentially estimated by counting gravitational wave cycles, so for a double neutron star coalescence 
it will be estimated to a precision that is typically b ette r than ~ 0.1% given that ~ 10 3 cycles are 
expected in the band of the detectors (see Berr y et al. 2Q14 for a recent study). Thus we can 
assume that M c h will be determined exactly, but we cannot necessarily assume that the signal will 
be strong enough to break the degeneracy and determine both masses independently. We can invert 
the equations to find the individual gravitational masses, given the unknown mass ratio q: 


M x = q 2 / 5 (l + g)!/ 5 M ch 
M 2 = 9 - 3/5 (l+9) 1/5 M ch 


Therefore, if we infer a given chirp mass, there will be some uncertainty in the individual gravita¬ 
tional masses and thus in the total baryonic rest mass in the merger. 

This uncertainty, however, is only a few hundredths of a solar mass, at least for masses similar 
to what we see in double neutron star systems in our Galaxy. This is fundamentally because 
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neutron star masses are limited from below; for precisely measured masses the current lower limit 
is 1.25 Mq, but for the sake of argument let us say that the true limit is 1.2 Mq. Therefore, 
for a given chirp mass, there is not much room to change the total baryonic mass. As a specific 
example, suppose we consider a merger between stars with gravitational masses M\ = 1.35 Mq 
and M 2 = 1.45 M 0 , in the context of the parameter combination for EOS1 that maximizes the 
gravitational mass of a nonrotating star for successful production of short gamma-ray bursts with 
these masses. The chirp mass for this combination is M c h = 1.218 Mq, and the total baryonic rest 
mass for the two stars combined and this EOS is Mb ary tot = 3.078 M 0 . If we take the extreme 
that one of the stars actually has a gravitational mass of Mi = 1.2 Mq, then to keep the same 
chirp mass it is necessary that M 2 = 1.64 Mq. The total baryonic rest mass for this combination 
is then My jarytot = 3.136 Mq. Thus even for this extreme case the difference is less than 0.06 Mq. 
If the lower limit to the gravitational mass of a neutron star is 1.25 Mq instead of 1.2 M 0 then the 
maximum total baryonic mass for this case drops to 3.107 Mq. If the chirp mass is smaller (as it 
will be if the gravitational masses of the neutron stars are lower) then the correction will be even 
less. We conclude that gravitational wave measurement of just the chirp mass from a coalescence 
coincident with a gamma-ray burst will place strong constraints on the total baryonic rest mass for 
a given EOS. 

If the maximum mass of a nonrotating neutron star is towards the high end of what we infer 
(say, 2.2 Mq) then mergers between lower-mass neutron stars will produce a remnant that will 
not collapse quickly unless angular momentum is actually removed from the system rather than 
redistributed. As we discussed, this likely requires the production of a very strong poloidal magnetic 
field within tens of milliseconds. If this does not happen in most cases, then many mergers could 
fail to produce short gamma-ray bursts; it would be interesting to know the observed properties 
of such an event. It would also suggest that estimates of the gravitational wave detection rate 
of double neutron star mergers based on the short gamma-ray burst event rate are conservative, 
because only some fraction of coalescences lead to bursts. 

Another possibility is that many short gamma-ray bursts are actually produced by the coales¬ 
cence of neutron stars with black holes rather than neutron stars with neutron stars. The larger 
chirp masses of such events means they will be detectable to greater distances than double neutron 
star mergers. Thus for a given observed gamma-ray burst rate per volume, the gravitational wave 
detection rate would be increased by a factor of several. 

It is important to note that, although in the context of our parametrizations an EOS that 
predicts a mass in excess of our maximum would not lead to short gamma-ray bursts, the converse 
is not necessarily true: a maximum nonrotating mass below our limit does not guarantee burst 
viability. For example, consider a 1.35 Mq — 1.45 Mq merger, and let us use EOS1 with a = 
13.3 MeV, a = 0.51, b = 4.1 MeV, j3 = 2.3, p t hresh,i = 1.85/o s , c s = 0.75c, and p t h re sh,2 = 3 p s , 
but (unlike in our standard parametrization) we set c s = c/2 above pthresh, 2 , then the maximum 
gravitational mass of a nonrotating star is M max = 2.06 Mq. This is well below our threshold for 
1.35 Mq — 1.45 Mq mergers with remnants whose spin is limited only by mass shedding. However, 
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the total baryonic rest mass of the stars, 3.075 Mq, is less than the maximum 3.108 Mq that could 
be supported by a uniformly rotating star, and therefore a short burst would not occur with this 
combination. Thus individual EO Ss sh ould b e tested ag ainst the short gamma-ray burst criterion 
using a code such as rns fStergioulas & Friedman 1995). 


Finally, we note tha t we find tighter constraints M max ~ 2 — 2.2 Mq than are reported in 
the study of Fryer et_ah ((.2015), who give limits of 2.3 — 2.4 Mq. We believe that the primary 
reason for the difference is that Fryer et al. (2015) are concerned with the question of whether any 
significant fraction of double neutron star mergers will produce short gamma-ray bursts. Thus they 
concentrate on the high end of the neutron star masses that emerge from their population synthesis 
calculations, which means that they find larger upper limits than we find by focusing on three mass 
categories that we know exist in the Galaxy. 

Again, we emphasize that when ground-based gravitational wave detectors see bursts with 
mergers, the lowest mass example of a successful burst will set the strongest limits on the maximum 
mass of a nonrotating neutron star. To that end it will be helpful to know whether there are any 
electromagnetic signatures of rapid collapse that can be identified clearly and seen from a broader 
range of directions than the burst itself, because this would increase the otherwise small fraction of 
mergers detected usin g gravitational waves that can be evaluated using our argument. For example, 
it has been proposed ((Metzger fc Fernandez! [2014]; iKasen et al.l 120151 ) that prompt formation of a 
black hole will lead to nearly isotropic red emission from disk winds within several days, whereas 
delayed formation of a hole will produce bluer emission within roughly a day. It is also conceivable 
that there is a signature in the gravitational wave emission itself of prompt collapse, although this 
is likely to be at frequencies > 2 kHz where ground-based detectors are not especially sensitive. 


In summary, we show that current models of short gamma-ray bursts involving double neutron 
star mergers imply a strong upper limit to the maximum mass of nonrotating neutron stars. If some 
of these bursts come from stars towards the low mass end of what we see in our Galaxy, then with 
the parametrized equations of state we have explored the limit could be close to M max ~ 2 Mq. 
If the bursts only come from neutron star binaries near the high end of our Galactic sample then 
the maximum could be M max ~ 2.2 Mq, but in that case short bursts are likely to occur only 
for high-mass neutron star binaries, which means that gravitational wave detection rates inferred 
from burst rates need to be increased. In either case, direct detection of gravitational waves from 
neutron star binaries along with coincident gamma-ray bursts will constrain strongly the properties 
of cold matter beyond nuclear density. 
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Fig. 1.— Upper limits on the maximum gravitational mass of a nonrotating star, based on the 
criterion that to make gamma-ray bursts, double neutron star mergers must produce an object 
that collapses when it becomes uniformly rotating. Red symbols refer to EOS1 and blue symbols 
refer to EOS2 (see text for descriptions of these EOSs). The three types of mergers we consider 
are of gravitational masses 1.25 M 0 and 1.35 M 0 ; of 1.35 M 0 and 1.35 M 0 ; and of 1.35 M 0 and 
1.45 Mq, to represent the range of masses in known double neutron star binaries. Filled squares 
show the limits obtained when we assume that the rotation of the merged remnant is limited only 
by mass shedding, which is our standard case and which is consistent with existing simulations. 
To explore the probably unrealistic situation in which there is very rapid angular momentum loss 
from the system, we also consider cases in which there are stricter angular momentum limits: 
crosses show the masses obtained when we assume that the more restrictive of mass shedding or 
T/W <0.1 limit the rotation of the remnant, and the open symbols are for nonrotating remnants. 
The dotted horizontal line shows the maximum currently known gravitational mass 2.01 M 0 for a 
neutron star. This figure demonstrates that if short ganmra-ray bursts are produced by the merger 
of double neutron star systems comparable to the ones we see in our Galaxy, then the upper mass 
limit for nonrotating neutron stars is constrained tightly. 




